Odd-frequency superconducting pairing in topological insulators 
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We discuss the appearance of odd-frequency spin-triplet s-wave superconductivity, first proposed by Berezin- 
skii [JETP 20, 287 (1974)], on the surface of a topological insulator proximity coupled to a conventional spin- 
singlet s-wave superconductor. Using both analytical and numerical methods we show that this disorder robust 
odd-frequency state is present whenever there is an in- surface gradient in the proximity induced gap, includ- 
ing superconductor-normal state (SN) junctions. The time-independent order parameter for the odd-frequency 
superconductor is proportional to the in- surface gap gradient. The induced odd-frequency component does not 
produce any low-energy states. 

PACS numbers: 74.45,+c, 74.20.Rp, 74.50.+r 



I. INTRODUCTION 

Topological insulators (TIs) are a new class of materials 1,2 
with an insulating bulk but with a conducting surface state. 
The surface state has its spin locked to the momentum in 
a Dirac-like energy spectrum. Superconducting TI surfaces 
have received a lot of attention recently, 3,4 since it was pre- 
dicted that Majorana modes appear in e.g. superconduct- 
ing vortex cores and at superconductor-ferromagnet (SF) 
interfaces. 5,6 Experimentally, both superconducting transport 7 
as well as the Josephson effect 8 have already been demon- 
strated in TIs proximity-coupled to conventional supercon- 
ductors. Despite this large interest, relatively little attention 
has been paid to the superconducting state itself. In addition 
to the standard proximity effect, one could expect the spin- 
orbit coupling in TIs to lead to significant modifications and 
produce novel superconducting states that are not easily ac- 
cessible in conventional superconductors. For example, it has 
already been demonstrated that an effective p-wave pairing is 
induced when the TI is proximity-coupled to a conventional 
s-wave superconductor.^*^ 

Quite generally, the superconducting pair amplitude, being 
the wave function of the Cooper pairs, needs to obey Fermi- 
Dirac statistics. This leads to the traditional classification into 
even-parity (s, d, ...) spin-singlet and odd-parity (p, f, ...) 
spin-triplet pairing. It has also been shown that the pair ampli- 
tude can be odd in frequency. 11,12 Odd-frequency spin-triplet 
s-wave pairing has been found in spin-singlet s-wave SF junc- 
tions due to spin-rotational symmetry breaking and it ex- 
plains the long-range proximity effect in these junctions. 13,14 
Very recently, the same magnetic field induced odd-frequency 
pairing has also been found in TIsJ^ In superconductor- 
normal metal (SN) junctions translational symmetry break- 
ing instead generates odd-frequency spin- singlet odd-parity 
component s 16 ! 17 . However, the odd-parity limits the odd- 
frequency pairing to ballistic junctions^ 

In this article we show with a simple analysis that the ef- 
fective spin-orbit coupling k • a on the TI surface immediately 
induces odd-frequency spin-triplet s-wave correlations, even 
in the absence of a magnetic field. The s-wave nature of the 



odd-frequency component makes it robust against disorder, in 
sharp contrast to normal SN junctions. The odd-frequency 
correlations appear whenever there is an in- surface gradient 
in the proximity-induced spin-singlet s-wave pairing, with the 
odd-frequency order parameter directly proportional to the in- 
surface gradient. We numerically calculate the odd-frequency 
response in several superconducting two-dimensional (2D) TI 
systems, including SN, SS' junctions, and in the presence of 
surface supercurrents. These results point to an important 
missing component in the discussion on the role of proximity 
induced superconductivity in TIs and the odd-frequency com- 
ponent ought to be included in the study of low energy states 
and Majorana fermions in TIs. We also discuss experimental 
consequences of the odd-frequency pairing. We find that the 
analytic form of the odd-frequency response does not re- 
sult in low energy states, which previously has been intimately 
linked to the appearance of odd-frequency components 
but the predicted gapped spectrum could allow detection of the 
odd-frequency component with local tunneling probes. The 
spin-triplet pairing will further produce a finite Knight shift 
in nuclear magnetic resonance (NMR) or muon spin-rotation 
measurements. 



II. ANALYTIC DERIVATION 

We start with an analytic calculation that illustrates the ap- 
pearance of an odd-frequency component. The Hamiltonian 
that describes the contact between a TI and a conventional s- 
wave spin-singlet superconductor (SC), see Fig Uta), can be 
written as H = H^i + i^sc + Ht'. 

k,a,/3 

^SC = e ( k ) d a,k d «,k + A W"/3 d a,»4,i + Ha 

k,a,/3 i,a,/3 

H T = J2T i clid a , i + n.c.. 
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H T i is the Hamiltonian describing the TI surface state at 
y < 0, Hsc is the Hamiltonian for the SC electronic states 
at y > 0, and Ht describes the tunneling between the SC and 
TI. We use the low energy approximation around the Dirac 
point for the TI dispersion. Hsc is of the standard form, 
where we explicitly allow for position dependence through 
the site index i, where then A = A(i) a p = e a f3A{i) is 
the matrix in spin space, which we choose to be spin-singlet. 
We assume that the kinetic energy e(k) in the SC is a sim- 
ple band and that the tunneling matrix element Ti is nonzero 
only for neighboring sites at the TI-SC interface. We further 
assume that A(%) is dependent on the in-surface x-coordinate 
and approximate it by a linear slope A(i) = A + (ia)^\o 
where a is the unit cell size. Superconducting correlations 
in the TI will be induced by pairs tunneling into the TI from 
the superconducting side. Therefore, any induced pair am- 
plitude in the TI will be proportional to A(i) and its deriva- 
tives at the interface. To derive the induced pair component 
on the TI surface, we evaluate the anomalous Green's func- 
tion F T i,a</3(k, k') = -z(T r Ca,(T,k)c/3(0,k')). Using stan- 
dard methods we find it to be proportional to 



FtiKIm) = -\T\ 2 Y,G\uj n \i,j)F(j,l)G°(uj n \l,i), 
hi 



(2) 



where G°(u;|k) = (k • a - iuj n )/(k 2 + cj£) is the Matsub- 
ara Green's function for the electrons in the TI, F(uu n , k) = 
A(i)[uj 2 + e 2 (k) + A(i) 2 ]~ l is the anomalous Green's func- 
tion for a conventional superconductor, where we assume 
slow variations of A(i). By going to momentum space and 
using the /c-space expansion A(k) = A £k,o + ^f 1 lo^k^, we 
can rewrite the nontrivial part of the induced pair amplitude 
on the surface of the TI as 



F TI (u n \i = 0)=J2 : 

k 



-z|T| 2 ^A| G Q K|k)a kcc G Q K|-k) 
2[a;2+ e 2( k ) + A2( )] 

ITiWa.Aio 



2[o;2+ e (k)2 + A2(0)](a;2+k2)2- 



(3) 



This result indicates that the Josephson tunneling into the TI 
will immediately induce an odd-frequency spin-triplet even- 
momentum superconducting component. 11 It is the character- 
istic spin-momentum locking in the Dirac surface spectrum 
that induces spin-triplet odd-frequency pairing in the pres- 
ence of translational symmetry breaking. The s-wave na- 
ture of this pairing guarantees robustness against disorder. 
The situation here is different from the normal metal case 
where only odd-frequency spin- singlet odd-parity correlations 
are induced by translational symmetry breaking ilSJl To eval- 
uate the on-site odd-frequency component we will assume 
that E F of the conventional supercoductor is the dominant 
scale and we ignore higher order terms in 1/Ep. The local 
on-site amplitude on the TI surface is given by an integral 
over the momenta and is proportional to F T i(uj n \i = 0) ~ 
\T\ 2 uj n a z d x A\ /(E^\uj n \ 2 ). Interestingly, 1/w dependence 



has also been reported in heavy fermion compounds. 21 The 
particular form of the gap function allows us to introduce an 
order parameter, i.e. the inherent parameter that is constant at 
equal time in the odd-frequency superconductor, see e.g. [22l 
The odd-frequency order parameter in our case is proportional 
to 



d T F Tl (r\i)\ 



\T\ 2 a z 



dA 



dx 



dA 

dx 



(4) 



This proportionality to the in-plane gradient of the s-wave gap 
can be tested and we indeed find that d T F^\ (r\i) |o is tracking 
the spatial gradient of the gap, see Figs.[TJf) and[2b). 

One of the observable consequences of odd-frequency 
superconductivity is usually the appearance of sub-gap 
states, 16,17,19,20,23 or even a low-energy continuum associated 
with a gapless nature of the state. 12,22 However, we find 
here that the particular structure of the odd-frequency gap 
~ l/|cj n | results in no intragap states at the lowest energies 
and thus this odd-frequency gap state is fully gapped. In- 
deed, after analytic continuation from the Matsubara axis, the 
local density of states (LDOS) is proportional to N(E) ~ 
Re[E/y/E 2 - A(E) 2 ] ~ Re[l - C 2 /£ 2 ]-\ which van- 
ishes at energies below the minigap induced by the tunneling 
C ~ \T\ 2 dA/dx. One has to keep in mind that these features 
will occur in addition to the LDOS features introduced by the 
induced conventional BCS pairing. 



III. NUMERICAL RESULTS 

The analytical results in Eqs. ([3]|4]) are derived for a 3D TI, 
but are equally valid for a 2D TI. To complement these results 
we explicitly calculate the odd-frequency spin-triplet super- 
conducting correlations in the Kane-Mele 2D TI:— 



KM 



(id) 



^c]ci + i\ i/ijc]a~c 



/ y I 



((hj)) 
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(5) 



where c\ is now the fermion creation operator on site i in the 
honeycomb lattice with the spin-index suppressed. Further- 
more, and denote nearest neighbors and next- 
nearest neighbors respectively, t is the nearest neighbor hop- 
ping amplitude, fi the chemical potential, A the spin-orbit cou- 
pling, and Vij = +1 (— 1) if the electron makes a left (right) 
turn to get to the second bond. We set the energy and length 
scales by fixing t = 1 and a = 1, respectively and choose 
A = 0.3, which gives a bulk band gap of 1. The influence 
of the SC can be described by an effective on-site attractive 
potential Ui acting at the TI edge when it is in proximity to a 



(6) 



We solve H = Hkm + #a self-consistently for the 
spin-singlet s-wave mean-field order parameter A(i) = 
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—Ui(ciiCi^). We further introduce the odd-frequency spin- 
triplet s-wave pairing correlations: 26 

F t °(T\i) = (c it (r) Cil (0) + c ii (r)c it (0)>/2, (7) 
F}(t\i) = (c iT (r)c it (0) - c^(r)c 4 (0))/2 (8) 

where r is the time coordinate (at zero temperature). Eqs. 
[8]) contain all space and time information of the odd-frequency 
spin-triplet response. As required by the Pauli principle, 
F t (r = 0) always vanishes for a self-consistent solution 
of A, whereas the time derivative at equal times d T F t \o 
defines the odd-frequency order parameter, in analogy with 
Eq. d3]). Also, since Hkm commutes with a z , the m = 1 
spin-triplet projection F^ is identically zero. We also de- 
fine a time-dependent quantity for s-wave spin-singlet pairing: 
F a [r\i) = ((c^(r)c it (0)-c it (r)c^(0)))/2, which is related 
to the order parameter through A$ = —UiF s (r = 0|i). 

A. SN junction 

From Eq. ([3]) we know that non-zero odd-frequency spin- 
triplet correlations require a gradient in the superconducting 
order parameter along the TI edge, i.e. ^ needs to be finite. 
This is the case e.g. at any step edge in a TI proximity coupled 
to a SC, but a more striking example is a SN junction along the 
TI edge, as schematically pictured in Fig.QJa). In the S region 
of the TI, we apply a constant U, and, since the SC provides an 
ample source of charge, we also set /is > Mn- In Fig. [Tib) we 




FIG. 1: (Color online) (a) SN junction along the edge of a 2D TL (b) 
Even-frequency pairing |F S (0)|, (c) odd-frequency pairing |F t °(l)|, 
and (d) |F t (6)| in a 2D TI SN junction with /x s = 0.3, /x N = 0, and 
U = 2.2, giving As = 0.16. (e) Maximum height Jif (black, left 
axis) and average width wf (red, right axis) as function of r for |F t ° | 
(crosses) and \d x F s \ (line), (f) Magnitude of odd- frequency order 
parameter d T Ft |o (crosses) and |A scaled by a factor of 0.4 (line) 
for A s = 0.16 (black) and A s = 0?30 (red). 

plot the magnitude of the self-consistent spin- singlet pairing 



amplitude F s (r = 0), which shows proximity -induced super- 
conducting pairing in the N region, as well as an inverse prox- 
imity effect (depletion of Cooper pairs) on the S side of the 
junction. The odd-frequency response is shown in Figs. QJc- 
e), with F t only non-zero in the SN interface region where the 
gradient of A is finite. For small r, F t is sharply peaked at 
the interface with the peak height h p rapidly increasing with 
r, as seen in Fig. Die). For r > 4 (small frequencies uo n ) 
the height is approximately constant, but the peak instead be- 
comes broader, with the average width wp, defined as the ra- 
tio of the total weight of the peak to its height hp, increasing 
roughly linearly with r. In Fig. [TJe) we also plot the height 
and width of the gradient of the spin-singlet response d x F s 
for a direct comparison. For large r, F t is directly related 
to d x F s , supporting the analytic result in Eq. ([3]), whereas at 
small r the direct relation breaks down for the peak height. In 
Fig - DID we finally compare the odd-frequency order param- 
eter d T F t |o with As predicted in Eq. (|4]), we find these 
two quantities to be directly proportional. 

The odd-frequency pairing in a SN TI junction is strik- 
ingly different from that of a conventional SN junction. The 
odd-frequency response in the TI case has even-parity spin- 
triplet symmetry, which is robust against impurity scattering 
and can thus survive even in the diffusive regime, 18 in con- 
trast to the odd-parity spin- singlet symmetry in a regular SN 
junction J^li Related to this disorder robustness, we also find 
that F t is insensitive to any Fermi level mismatch at the SN 
interface, created by using different chemical potentials in 
N and S. In a normal SN junction, the odd-frequency pair- 
ing quickly disappears when the transparency of the interface 
is reduced. 17 Finally, to investigate the influence of the odd- 
frequency pairing on the low-energy spectrum, we study SS' 
junctions, where the two S regions have different pair poten- 
tials U. We find no evidence for a reduced gap or intragap 
states in both sharp and extended interface junctions. 27 This 
expands our analytical low-energy result to exclude any intra- 
gap states, and is again in sharp contrast to conventional SN 
and SF junctionsi 16 i 17 i 19 i 20 i 23 



B. Supercurrent 

Odd-frequency correlations can also be generated in a ho- 
mogenous TI-SC system if a supercurrent I is applied in- 
surface, since the current is proportional to the gradient of 
the phase of the superconducting order parameter. We model 
such a system by setting A = \A\e zkx , with k being the 
(fixed) phase winding which is proportional to the current, 
and solve self-consistently for A. We find that F s (r) = 
C 3 (r)e< kx+6a ^ and F t = C t {r)e^ kx+et ^\ i.e. both F s 
and F t have the same phase winding k for all r, but the ampli- 
tudes C and phase off-sets 6 are dependent on r. In Fig. [2a) 
we plot the (position-independent) magnitudes of F t and d x F s 
as function of r for two different values of k. Similar to the SN 
junction, F t increases rapidly for small r. For small currents 
F t continues to increase even for large r, even though less 
steeply. This is in contrast to d x F s which, apart from small os- 
cillations, decreases with increasing time. For larger currents 
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FIG. 2: (Color online) (a) Magnitude (height) of the odd-frequency 
pairing F t ° (crosses) and d x F s (line) as function of r for phase wind- 
ings k = 7r/50 (black) and 7r/10 (red) along the TI edge. Inset shows 
schematically a supercurrent / along the edge of a 2D TI. (b) Mag- 
nitude of the odd-frequency order parameter <9 T F t °|o (crosses) and 
(line) as function of k. Inset shows their ratio. Here U = 2.5, 
/x = 0.3, giving A s = 0.30. 

we see a sharp downturn in F t at intermediate time values. Be- 
fore this downturn we find that F t is directly proportional to k, 
as expected from its relation to d x F s in Eq. ©. However, for 
times past the downturn, F t instead decreases with increasing 
k. The phase off- set parameters S and 6 t also tracks each 
other (with the expected tt/2 shift) before the downturn in F t , 
but past this time the correlation between them is lost. Thus, 
we find that F t and d x F s are only tracking each other in a fi- 
nite r-range, which set by the current. In Fig. Ob) we focus 
on the behavior at r = 0. Both the odd-frequency order pa- 
rameter d T F®\o and ^ are linear in k, with their ratio being 
a constant for all currents. This shows that the odd-frequency 
order parameter is always directly proportional to the current 
(phase winding), in agreement with Eq. ©, and thus applying 
a supercurrent offers an experimentally easy way to tune the 
strength of the odd-frequency pairing. 
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FIG. 3: (Color online) (a) Odd-frequency pairing 1^(1)1 in the SN 
junction in Fig. Q] with Xr = 2A. (b) \Ft | along the TI edge and (c) 
at the SN interface into the TI for r = 0.5 (dashed black), 2 (black), 
8 (dashed red), and 16 (red), (d) Spin components |F t +1 (l)| and (e) 
|F t - 1 (l)|ofF t 1 (l) = (F+ 1 -F t - 1 )/2. 



F^ 1 is generated whenever ^ is non-zero, i.e. for a finite 
order parameter gradient perpendicular to the surface of the 
TI. Since the low-energy density of states varies dramatically 
between the surface and the bulk of a TI, there will always 
be a strong such gradient for proximity-induced superconduc- 
tivity in a TI. In the bulk of the S region F^ 1 = Ff, and 
thus the criterion for F} to be non-zero is the same as for FP, 
i.e. a finite gradient in-surface gradient In terms of the d- 
vector for the d T F f | r =o odd-frequency pairing, we find that it 
is always real, yielding a unitary spin-triplet state. 



C. Rashba spin-orbit coupling 

So far we have only discussed m = spin-triplet pair- 
ing. Equal spin-triplet pairing (m = 1) is produced at the 
TI edge if a term misaligned with the spin-quantization axis 
is present in the Hamiltonian. This happens for Rashba spin- 
orbit coupling—: 

H R = iX R 4( s x <Uj)zCi, (9) 

which is present when z — » — z symmetry is broken. Here 
dij is the unit vector along the bond between sites j and i. 
Hkm + Hr is still in the non-trivial topological phase with 
a single Dirac cone at the edge for Xr < 2\^3X. We find 
that F t ° decreases only slightly when introducing a finite Xr, 
and even for = 3 A, we have F® > F} . For a SN junc- 
tion we find that F^ is localized to the interface region, hav- 
ing a similar x dependence as F t °, see Figs. Oa-c). Here we 
also see that F} extends somewhat farther into the TI, espe- 
cially for larger r. We further analyze F} by studying its two 
spin parts F^ 1 = (Q ) ± cr (r)Q 5 ± cr (0)). Surprisingly, F^ 1 is 
not only non-zero at the SN interface, but in the whole S re- 
gion, as shown in Fig.[3td-e). In fact, we find that a non-zero 



IV. SUMMARY 



In summary, we have shown that odd-frequency spin-triplet 
5-wave pairing is present in a TI proximity-coupled to a con- 
ventional 5-wave spin-singlet superconductor in zero mag- 
netic field. The time-independent order parameter for the odd- 
frequency pairing is proportional to the in-plane gradient of 
the induced s- wave gap. This disorder robust odd-frequency 
response is an immediate consequence of the spin-momentum 
locking in the TI surface state. We have explicitly demon- 
strated the occurrence of odd-frequency correlations not only 
in SN and SS' junctions at a 2D TI edge, but also when a su- 
percurrent is applied along the edge. In terms of experimen- 
tal observables, we find no evidence of subgap states in the 
presence of odd-frequency pairing, due to its particular fre- 
quency dependence. The gapped LDOS could allow the de- 
tection of the odd-frequency component with local tunneling 
probes. Furthermore, the spin-triplet component could pro- 
duce a Knight shift in NMR or muon spin-rotation measure- 
ments. 
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V. SUPPLEMENTARY MATERIAL 

In this supplementary material we provide numerical data 
showing the absence of subgap states in the presence of odd- 
frequency pairing in a topological insulator (TI) induced by 
a gradient in the proximity-induced conventional s-wave su- 
perconducting state. In order to do so accurately we explic- 
itly model the microscopic interface between a superconduc- 
tor (SC) and a two-dimensional (2D) TI displayed in Fig. 01 
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FIG. 4: (Color online) Microscopic details of the TI-SC interface 
with TI sites a (dark) and SC sites bi (light). 

The total Hamiltonian is H = Hkm + Hsc + Hf, where 
Hkm defines the TI and is given by Eq. (5) in the main text, 
whereas 

H SC = -f £ b\ a b ia -UJ2 bfablbn (10) 
H~ t = -i J2 4^+H.c, (ID 

defines the SC and the coupling between the TI and the SC, re- 
spectively. The SC is defined on a square lattice with nearest 
neighbor hopping t and an on-site spin-singlet s-wave pair- 
ing from an attractive Hubbard U term. The coupling be- 
tween the TI and the SC is by a tunneling element t acting be- 
tween nearest neighbors across the interface. We treat Eq. (TTOb 
self-consistently within mean-field theory by using the self- 
consistency condition A (i) = —U(biibrf). 

In Fig. [5] we show a typical interface when the pairing 
potential U is set to vary along the interface, i.e. in the re- 
direction, in order to produce a gradient along the TI edge, 
but is constant along the ^-direction. Region A has a constant 
U = 2.5 giving A = 0.60 in the bulk of the SC and region 
C has a similarly constant U = 4 giving A = 1.38. In the 
intervening region B there is a linear rise of U between these 
two values. Finally, the interface is put on a cylinder, making 
a sharp S(U = 2.5)-S(J7 = 4) interface at D. In Fig. 0a) 
we clearly see how the magnitude of the induced s-wave pair- 
ing Fjj in the TI reflects the change in U along the interface. 
Figure [3b) displays the magnitude of the odd-frequency pair- 
ing order parameter <9 r F t °| r= o, which is only non-zero in the 



B and D regions. The odd-frequency pairing leaking into the 
SC is at least an order of magnitude smaller and we can not 
deduce any physical consequences in the SC from this back 
action. Figures |3c)-(f) shows the corresponding local den- 
sity of states (LDOS) in the A-D regions. Starting from the 
far left, these LDOS plots show the unperturbed left TI edge, 
with a constant DOS due to the one-dimensional surface Dirac 
cone. When we see the constant TI bulk gap of 1 before the 
right TI edge which is gapped by the induced superconduc- 
tivity from the SC. The gap in the bulk of the superconductor 
equals A which varies significantly in x, but the induced gap 
in the TI surface state varies much less. Most notably, we see 
no evidence of any subgap states in the C and D region where 
odd-frequency pairing is present. In fact, the right TI edge 
in the C and D regions looks remarkably similar to a simple 
interpolation between the edge states in regions A and B. 
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FIG. 5: (Color online) TI-SC interface with U — 2.5 in the A region 
(18 sites), U = 4 in the C region (18 sites), linearly varying between 
these values in region C (16 sites) and an atomically sharp drop/rise 
in region D. (a) Magnitude of even-frequency pairing F s and (b) 
odd- frequency pairing d T Ft\ T =o in the TI. (c)-(f) LDOS plots with 
grayscale limits 1 (black) and (white) states/site/energy in regions 
A-D. Dotted line mark the Fermi level. Here /x = 0, A = 0.3, and 
t = 0.9. 



